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Abstract 

In this paper we consider a model of quantum computation based on n atoms 
of laser-cooled and trapped linearly in a cavity and realize it as the n atoms Tavis- 
Cummings Hamiltonian interacting with n external (laser) fields. 

We solve the Schrodinger equation of the model in the case of n=2 and construct 
the controlled NOT gate by making use of a resonance condition and rotating wave 
approximation associated to it. Our method is not heuristic but completely math- 
ematical, and the significant feature is a consistent use of Rabi oscillations. 

We also present an idea of the construction of three controlled NOT gates in the 
case of n=3 which gives the controlled-controlled NOT gate. 
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1 Introduction 



Quantum Computation (or Computer) is a challenging task in this century for not only 
physicists but also mathematicians. Quantum Computation is in a usual understanding 
based on qubits which are based on two level systems (two energy levels or fundamental 
spins) of atoms, See jlj as for general theory of two level systems. 

In a realistic image of Quantum Computer we need at least one hundred atoms. How- 
ever, then we may meet a very severe problem called Decoherence which destroy a su- 
perposition of quantum states in the process of unitary evolution of our system. At the 
present it is not easy to control the decoherence. See for example |2I as an introduction. 

An optical system like Cavity QED may have some advantage on this problem, there- 
fore we consider a quantum computation based on Cavity QED. As an approximate 
model we realize it as the n atoms Tavis-Cummings Hamiltonian interacting with n ex- 
ternal (laser) fields. As to the Tavis-Cummings model see P]. To perform the quantum 
computation we must first of all show that our system is universal j3] . To show it we must 
construct the controlled NOT operator (gate) explicitly in the case of n = 2, [Sj. 

For that we must embed a system of two-qubits in a space of wave functions of the 
model and solve the Schrodinger equation. In a reduced system we can construct the 
controlled NOT by use of some resonance condition and the rotating wave approximation 
associated to it. Then we need to assume that the coupling constants are small enough 
(the weak coupling regime in the title). 

Next we want to construct the controUed-controUed NOT operator in the case of 

n = 3. For that purpose the construction of three controlled NOT gates is required ^ 

because three atoms are trapped linearly in the cavity, and we present an idea toward 

explicit construction. If this point will be overcome our system of quantum computation 

may become complete. 

^In the study of Cavity QED Quantum Computation this (important) point is missed 
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2 A Model Based on Cavity QED 

We consider a quantum computation model based on n atoms of laser-cooled and trapped 
linearly in a cavity and realize it as the n atoms Tavis-Cummings Hamiltonian interacting 
with n external (laser) fields. This is of course an approximate theory. In a more real- 
istic model we must add other dynamical variables such as positions of atoms and their 
momenta etc. However, since such a model is almost impossible to solve we consider a 
simple one. 

Then the Hamiltonian is given by 

n 

J2 hj (af )e^(^^*+<^^) + aj-)e-*(^^*+'^^)) ® 1 (1) 

where u is the frequency of radiation field, A the energy difference of two level atoms, a 
and are annihilation and creation operators of the field, and g a coupling constant, Qj 
the frequencies of external fields which are treated as classical fields, hj coupling constants. 



and L = 2". Here crj , cr] and a] are given as 



Lcxc ^.(+) ^(-) ^^A ^(3) 



a)' = 12® ■■■(g) h^cTs^h® ■■ -(8)12 U - position) e M(L, C) (2) 
where s is +, — and 3 respectively and 



01 /oo fio /lo, 

^- = > ^3 = > l2 = ■ (3) 



In the case of n = 2 (which is the target through this paper) see the figure 1. Here 
we state our scenario of quantum computation. Each external field generates a unitary 
element of the corresponding qubit (atom) like a®)b where a, b G f/(2), while an photon 
inserted generates an entanglement among such elements like J2j O'j ® bj. As a whole we 
obtain any element in U{4). 

Here let us rewrite the Hamiltonian (^. If we set 

^.=E-r^ ^-=E-r^ ss=w-?\ (4) 

j=i i=i ^ i=i 
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Figure 1: The dotted line means a single photon inserted in the cavity and two curves 
mean external (laser) fields (which are treated as classical ones) subjected to atoms 



then can be written as 

H = cull's) a^a + AS3 1 + g (S+ S a + S- S a^) + 



J2 hj (af ^e'(^^*+^^) + aj"^e-'(^^*+'^^)) S 1 = Hq + V{t), (5) 
i=i 

which is relatively clear. Hq is the Tavis-Cummings Hamiltonian and we treat it as an 
unperturved one. We note that {5+, S-, S3} satisfy the su(2)-relation 

[^3, S^] = S^, [^3, S.] = -S., [S^, S.] = 2S3. (6) 

However, the representation p defined by 

p{a+) = S+, p(o-_) = 5'_, ^(0-3/2) = 5*3 

is a full representation of su{2), which is of course not irreducible. 
We would like to solve the Schrodinger equation 

ijU = HU= {Ho + V) U, (7) 

where [/ is a unitary operator. We can solve this equation by using the method of 
constant variation. The equation i-^U = HqU is solved to be 



where N = a^a is the number operator and Uq a constant unitary. Here we have used the 
resonance condition 

a; = A (8) 



, see for example 0. By changing Uq i — ^ Uo(t) and substituting into ((Tj) we have the 
equation 

dt 

after some algebras. We would like to calculate the right hand side of © explicitly, which 
is however a very hard task due to the term Q-'^^9{s+(^a+S-m'^) ^ j^. ^^^^ been done only for 

n = 1, 2 and 3 as far as we know, [Zj, [0]. The case n = 1 which is just the Jaynes- 
Cummings model is not interesting from the point of view of a quantum computation, so 
we restrict to the case n = 2 in the following. 

First let us write down each term (jHl). From the result in [Hj and some algebras we 
have 

^ |f±|/(Ar + 1) + 1 -ik{N+l)a -zk{N+l)a ^f{N + l)a^ ^ 



-ik{N)a^ f{N) + l f{N) -ik{N)a 

-ik{N)a^ f{N) f{N) + l -ik{N)a 

\^J{N-l)a)^ -ik{N-l)a^ -ik{N - l)a) J^/(iV - 1) + 1 j 



(10) 



where 



and 



f{N) 



-1 + cos (t^^2(2iV+ 1) 



k{N) 



sin (tgsj2{2N + 1] 



2{2N +1] 



V{t) U-^^^-^ ® e-**"^^ 



}i^Qi{{^2+'-')t^-4>2} I^^Qi{(Sli+ui)t+<t,i} Q 

}l^Q-i{{^2+i^)t+tt>2} Q Q ^^gi{(f2i+i^)i+</>i} 

l^^Q-i{{^l+i^)t+tt>l} Q Q ^^gj{(C2+w)t+</'2} 

/^^g-i{(f^i+c^)i+0i} j^^^--i{(Q.2+uj)t+4>2} Q 



Therefore we can calculate the term 



® 1. 
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from ()10|) and (fTT|) . However, we omit the explicit form because of being too complicated. 

Next let us go to a quantum computation based on two atoms of laser-cooled and 
trapped linearly in a cavity. 



3 Quantum Computation 



Let us make a short review of two-qubits. Each element can be written as 



^ = a++|+) ® 1+) + a+_|+) ® |-) + a-+|-) ® |+) + a__|-) ® |-) 



with two bases |+) and |— ) and |a++p + |a+_p + |a_4-p + |a_ 



1. Here if we identify 



1+) 



1 





1 



then ip above becomes 



(a ^ 



a |_ 

a 



(13) 



How do we embed two-qubits in our quantized system ? It is not known at the 
moment, which will depend on some method of experimentalists. Therefore let us consider 
the simplest one like 



m)) 



®|0), 



(14) 



a+_(t) 
a_+(t) 

where |0) is the ground state of the radiation field (a|0) = 0). We note that in full theory 
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we must consider the following superpositions 



\m) = E 



n=0 



® \n) 



a-+,n(t) 

V a__,„(t) y 

as a wave function, which is however too complicated to solve. 

To determine a dynamics that the coefficients 
stitute (^U into the equation 



will satisfy we sub- 



(15) 



See Appendix for the full calculations. The equation is not satisfied under the restrictive 
ansatz (fT^ . However, excited states |1), |2), |3) which have no corresponding kinetic 
terms contain the coupling constants hx and /i2, so the equation is approximately satisfied 
if they are small enough (namely, in the weak coupling regime in the title). 
Therefore the (full) equation is reduced to the equations of } at the 

ground state. 

. d 

f^^^r{in,+.)t+M |;(o) + ^/(0)/(l) + M0)M1)} + 
^^^,{^n,+.)t+M |i + ;(o) + ^/(i) + ^/(0)/(l) + A;(0)A:(1) > + 



h^e^^'^''^^-^'^^^^ \ 1 + /(O) + ^/(l) + ^/(0)/(l) + mk{l) \ + 



^ ^(Q) ^ ^/(0)/(l) + mk{l) 



(16) 



'dt''^- = 

/^ie-^«''^+")*+*^^ {/(O) + ^/(0)/(l) + A;(0)A;(1)} + 

1^ ^ ^(Q) ^ 2^(^) ^ ^/(0)/(l) + mk{l) 
hxeimi+-^)i+M |i + j(o)} + /i2e^«^^+^)*+^2>/(0) 



(17) 



d 

i—a- 
dt 



{l + /(O) + ^/(l) + ^/(0)/(l) + A;(0)fc(l)} + 



^^g-^{(n.+a;)*+<^.} 1^(0) ^ -/(0)/(l) + MO)Mi; 



hie 



i{(ni+uj)t+<f>i} 



/(0) + /i2e 



i{(Q2+t^)t+</'2} 



{l + /(0)} 



a++ + 
a , 



(18) 



I— a = 

dt 



h^e-^m^+'-)t+'Pi} {1 + /(o)} + /i2e-*^(^^+")*+^2^/(0)] a+_ + 



or in a matrix form 



d 
I— 
dt 



a+4t) 
a_+(t) 



a_+(i) 



(19) 



(20) 



where (i is the corresponding matrix element from the above equations. 

We obtained the system of complete equations, which is still comphcated. How do 
we solve it ? We use some resonance condition and the rotating wave approximation 
associated to it. Since 

/(O) = {-1 + cos (tgV2)}/2, /(I) = {-1 + cos {tgVe) } /2, 
k{0) = sin (tgV2) /V2, k{l) = sin (igVo) /VQ, 

the products /(0)/(l) and k{0)k{l) contain the term Q-^t9{V2+V6) |-,y ^j^g Euler formulas 
cos(^) = (e^^ + e-^^)/2, sin(^) = (e*^ - e-^^)/2i Noting 

we set a new resonance condition 



Qi+u;-{V2 + V6)g^0. 



:2i] 



All terms in except for the constant one Q'-{i^T-+'^-i^+^)9)t+'t>i} = q^<I>i contain ones like 
gi(te+a) _^ neglect all such oscillating terms (a rotating wave approximation). 

Then pn|l reduces to a very simple matrix equation 

/ 



a+_(t) 
a_+(t) 

The solution is easily obtained to be 



d 
I— 
dt 



24 















-«<P1 



^ a++(t) ^ 
a+_(t) 
a_+(t) 



exp 



t(y3- 1)M 
24 







V 











/ J 



cos(Q;t) ^^^sin(a)f:) 
-^sm(at) 

iQ-i<l>l . I ,\ 

^^sm(at) 



ie*"^! 



l+COS(at) 
2 



^sin(Q;t) 




l+COS(at) 
2 



-l+COS(at) 
2 

l+COS(at) 
2 



V 





( a++(0) ^ 







a++(t) 
a+_(t) 
a_+(t) 



( a++(0) ^ 
a+_(0) 
a_+(0) 

V «-(o) y 

^ ^ a++(0) ^ 
a+_(0) 
a_+(0) 



(22) 



(23) 



a+_(0) 
a_+(0) 

V «-(o) y 

where we have set a = ^^^ /ii. That is, we obtained the unitary operator U{t). In 
particular, if we choose to satisfying cos(ato) = —1 (sin(ato) = 0), then 



1 

1 



(24) 



-1 



9 



At this stage we use a very skillful method ^. That is, we exchange two atoms (jH]) in 
the cavity 

Exchange 



which introduces the exchange (swap) operator 

/ 



P 



1 



1 

1 



Multiplying ?7(to) by P gives 



PU{to) 



V 



(25) 



(26) 



This is just the controlled operator except for the overall constant —1 (an overall 
constant can be always neglected). From this it is easy to construct the controlled NOT 
operator, namely 

/ 



Cnot = (i2®vr)a.(i2®w^) 



1 



1 



v 



1 

1 



'^U{tQ) is imprimitivc in the sense of [H], so the main theorem in it says that our system is universal 
(namely, we can construct any element in C/(4)) . However, how to construct a unitary element explicitly 
is not given in |S] 
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where W is the Walsh-Hadamard operator given by 




(27) 



See for example As to a construction of W by making use of Rabi oscillations see jTHj. 
Therefore our system is universal jl], [Hj- 
A comment is in order. 

(a) In the equation (|20|l we can set another resonance condition in place of (|2T|l and obtain 
a unitary operator like U{t) in ()2H|1 . 

(b) In place of the ansatz ()14|) we can set for example 

/ ^ 

a+_(t) 





^ a++(t) ^ 



\m) 



\ 



® |0) 






a_+(t) 

V J 

Then we can trace the same line shown in this section and obtain a unitary operator 
under some resonance condition like (pT|l . This is a good exercise, so we leave it to the 
readers. 



4 Controlled-Controlled NOT Gate 

Our quantum computation model is based on n atoms of laser-cooled and trapped lin- 
early in a cavity, so we have another problem on the controlled NOT operators (of three 
types) when n = 3. 

Problem : Let us consider the case of three atoms in a cavity. How can we construct 
C-NOT (or C-unitary) operators for any two atoms among them ? 

See the figure 2. These constructions are very crucial in realizing quantum logic gates, 
for example, the controlled-controlled NOT gate shown as a picture 
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cc 



NOT 



or in a matrix form 



v 



1 

1 



The (usual) construction by making use of controlled NOT or controlled U gates is shown 
as a picture ([H], 
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-® — &■ 



where is a unitary matrix given by 



-4 



1 + i 1-i 
\^ 1 - i 1+i J 



1 

1 



CTi. 



However, we have not seen "reahstic" constructions in any references, so we must give 
the exphcit construction. 

To solve this let us state our idea. First we consider the construction of controlled 
NOT operator between the first and second atoms, namely 

C-NOT 



X)— = a 



NOT 



Our strategy is as follows. 

(i) We move the third atom from the cavity. 

(ii) We insert a photon in the cavity as two atoms interact with it and subject laser fields 
to the atoms, and next exchange the two atoms, which gives the controlled NOT operator 
as shown in the preceding section. 

(iii) We return the third atom (outside the cavity) to the former position. 
See the figure 3. 

If an influence of the "getting the third atom in and out" on the states space is small 
enough (namely, the unitary operator induced is near to the identity I4), then we certainly 
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Figure 3: The process to construct the controlled NOT gate between the first atom and 
second one for the three atoms in the cavity 
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Figure 4: The general setting for a quantum computation based on Cavity QED. The 
dotted hne means a single photon inserted in the cavity and all curves mean external 
(laser) fields (which are treated as classical ones) subjected to atoms 

obtain the controlled NOT gate (namely, Cnot ® I2) that we are looking for. Similarly 
we can obtain the remaining two ones. 

It is easy to generalize our idea to the n-atoms case. To perform a quantum com- 
putation we need to construct (many) controUed-controlled NOT gates or controlled- 
controlled unitary ones for three atoms among n-atoms, see §7. The method is 
almost same, so we leave it to the readers. See the figure 4. 

In principle, we can construct general quantum networks ^ 

By the way, a quick construction of controUed-controlled NOT gates is essential in 
general quantum networks |3]. 

We have given the exact form of evolution operator for the three atoms Tavis-Cummings 
model [Zj, therefore we can in principle track the same line shown in this paper and it may 
be possible to get the controUed-controlled NOT or many unitary gates directly (without 
combining many elementary gates like the construction of controUed-controlled NOT gate 
above). 

However, such a calculation for the three atoms case becomes very difficult (because 
we must treat 8x8 matrices at each step of calculations). We will attempt it in the near 
future. 

■^We must estimate an influence of the "getting atoms (which are not our target) in and out" on the 
whole states space, which is however difficult in our model. For that we must add in further terms 
necessary to calculate it 
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5 Discussion 



In this paper we constructed the controlled NOT operator in the quantum computation 
based on Cavity QED which showed that our system is universal. We also constructed 
the controlled-controlled NOT operator (under some assumption). Therefore we can in 
principle perform a quantum computation. 

We expect strongly that some experimentalists will check whether our method 
works good or not. 

See ^2] and their references for some experiments on Cavity QED (which may be related 
to our method). 

We conclude this paper by making a comment (which is important at least to us). 
The Tavis-Cummings model is based on (only) two energy levels of atoms. However, an 
atom has in general infinitely many energy levels, so it is natural to use this possibility. 
We are also studying a quantum computation based on multi-level systems of atoms (a 
qudit theory) ^1]. Therefore we would like to extend the Tavis-Cummings model based 
on two-levels to a model based on multi-levels. This is a very challenging task. 

Acknowledgment. We wish to thank Shin'ichi Nojiri for his helpful comments and sugges- 
tions, and thank the referees for careful readings and useful suggestions. K. Fujii wish to 
thank Gilles Nogues for teaching him some experimental facts on Cavity QED. 

Appendix 

One Qubit Operators by Classical Fields 

Let us make a brief review of theory without the radiation field, whose states space is 
only tensor product of two level systems of each atom. See the figure 5. 
The Hamiltonian in this case is 

H = E + (aj+)e*(^^*+<^^) + aj-)e-^(^^*+'^^)) | (28) 
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Figure 5: The n atoms in the cavity without a photon (in Figure 4) 



from where we have omitted the unit operator 1. Then 



/ 



E 



A 
2 



e 2 



e 2 





A 




2 




) 


2 




A 

2 






A 
2 



2 

"-J 2 



e 2 



e* 2 



/J 



(29) 



where 



e* 2 



V 



e 2 



and = I2 ® ■ ■ ■ ® I2 ® M (g) Is 



The wave function defined by i4:\^) = H\^) with pH|l can be written as a tensor product 



\^) = \iJi)®---^\iJn), (30) 

so if we define 



then it is easy to see 



|^) = (?7i®---®f/„)t|^), 



Th solution is easy to obtain 



i^(t)) = 0exp|-.ti \ )\^m- 
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Therefore, the solution that we are looking for is 



n ( 

= 

3=1 \ 

Last we note that 



e* 2 



exp I —it 



where 



exp < —it 



2 
h 



Xii Xi2 



Xu = COS I ^ + 



I sin t 



— i- 



02 



612 



Xi2 — X21 



X22 = COS 



-ih- 



sin ( t^/'i + 



02 



tf-^ + h^ 



sin (t./^ + /i2 



2 -/f + Z^^ 



1^(0)). (31) 



We can always construct unitary operators in f/(2) at each atoms by using Rabi 
osillations, see for example [TUj . 



Explicit Form of the Equation (II 5p 

Let us give the explicit form to (fTK|l for avoiding errors in the calculations. The left 
hand side of (fT3j) is 



z|a_+(t) 



® |0) 



(32) 



, while each of the right hand side becomes 



1-component = /i^e^«^^+'^)*+'^i> 



X 
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ia++ + -Ml)/(2) - 3/(1)^2)} |1) + «+- {/(O) + -/(0)/(l) + ^0)^1)} |0> 

+a_+ {1 + /(O) + ^/(l) + ^/(0)/(l) + fc(0)fc(l)} |0)] + 

-za++ + ^Ml)/(2) - ^/(l)fc(2)} |1) 

+a+_ |l + /(O) + + ^/(0)/(l) + Mo)Mi)} |o) 

(33) 



+a_+{/(0) + ^/(0)/(l) + MO)Ml) 

2-component = hie'^^^^+'^^^+'l'''^ x 
%/2a++ {^/(l) + k{l)k{2) + ^/(1)/(2)| |2) - za+_ {fc(0) - /(O)fc(l) + fc(0)/(l)} |1) 
-ia_+ {^(0) - ^(1) - /(0)A;(1) + fc(0)/(l)} |1) + a__ {1 + /(O)} |0)] + 

^2gj{(n2+a;)t+<;!.2} ^ 

V2a++ |a;(1)A;(2) + ^/(l)/(2)} |2) + za+_ + f{0)k{l) - k{0)f{l)} |1) 
+m_+ {/(0)A;(1) - M0)/(1)} |1) + a__/(0)|0)] + 
^^^-^{in,+.)t+Ma^^ {/(O) + ^/(0)/(l) + ^(0)^(1)} |0) + 

^^^-ii^n,+.)t+Ma^^ {1 + /(O) + ^/(l) + ^/(0)/(l) + M0)A:(1)} |0), (34) 

3-component = hie'^i^'+^)t+'Pi} x 
v^a++ {fc(l)A;(2) + ^/(1)/(2)| |2) + ia+_ {/(0)A;(1) - A;(0)/(1)} |1) 



3 

+ta.+ {k{l) + /(O)fc(l) - fc(0)/(l)} |1) + a-/(0)|0)] + 

x/2a++ {^/(l) + fc(l)fc(2) + ^/(1)/(2)| |2) - za+_ {fc(0) - k{l) - /(O)fc(l) + k{0)f{l)} |1) 
{A;(0) - /(0)A;(1) + ^(0)/(l)} |1) + a__ {1 + /(O)} |0)] + 
j^^^-ii^n,+^)t+M^^^ {1 + /(O) + ^/(l) + ^/(0)/(l) + A;(0)M1)} |0) + 
;,^e-H(n.+.)*+<^.}^^^ {/(O) + ^/(0)/(l) + M0)M1)} |0), (35) 

4-component = hie'^^^^+'^^^+'f'''^ x 
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-2x/6za++ {-fc(l)/(2) - 3/(1)^(2)} |3) + v^a+_ {-/(0)/(l) + ^(0)^(1)} |2) 
+v^a_+ 1^/(1) + ^/(0)/(l) + A;(0)M1)} |2) + m__MO)|l)] + 

-2v^za++ {^Ml)/(2) - ^/(1)M2)} |3) + v^a+_ + ^/(0)/(l) + ^0)^1)} |2) 



+v^a_+ {^/(0)/(l) + fc(0)A;(l) \ |2) + m__A;(0)|l) 



+ 



^^g-i{(Oi+i^)t+</.i} ^ 

+ ^fc(0)/(l) - 2/(0)fc(l)} |1) + a+- {1 + /(O)} |0) + a_+/(0)|0)] + 

' (36) 



m++ |fc(0) - + ^M0)/(1) - 2/(0)fc(l)} |1) + a+-/(0)|0) + a_+ {1 + /(O)} |0) 
after a long calculation by making use of ()12|1 . 
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